BAR4R-PUB-01/21 
SLAC-PUB-9012 



> 

(N 



O 



Study of CP-violating asymmetries in — ^ tt+tt , K~^7r decays 

B. Aubert,^ D. Boutigny,i J.-M. Gaillard/ A. Hicheur/ Y. Karyotakis/ J. P. Lees/ P. Robbe,^ V. Tisserand,^ 

A. Palano,2 A. Pompili,^ G. P. Chen,^ J. C. Chen,^ N. D. Qi,^ G. Rong,^ P. Wang,^ Y. S. Zhu,^ G. Eigen 

B. Stugu,"* G. S. Abrams,^ A. W. Borgland,^ A. B. Breon,^ D. N. Brown,^ J. Button-Shafer,^ R. N. Cahn,^ 

A. R. Clark,5 M. S. Gill,^ A. V. Gritsan,^ Y. Groysman,^ R. G. Jacobsen,^ R. W. Kadel,^ J. Kadyk,^ L. T. Kerth,^ 
Yu. G. Kolomensky,^ J. F. Kral,^ C. LeClerc,^ M. E. Levi,^ G. Lynch,^ P. J. Oddone,^ A. Perazzo,^ M. Pripstein,^ 
■ • ■ N. A. Roe,'^ A. Romosan,^ M. T. Ronan,^ V. G. Shelkov,^ A. V. Telnov,^ W. A. Wenzel,^ P. G. Bright-Thomas,*^ 
O ' T. J. Harrison,^ C. M. Hawkes,^ D. J. Knowles,^ S. W. O'Neale,^ R. C. Penny,^ A. T. Watson,^ N. K. Watson,^ 

T. Deppermann/ K. Goetzen/ H. Koch/ M. Kunze/ B. Lewandowski/ K. Peters/ H. Schmuecker/ M. Steinke/ 
^ ■ J. C. Andress/ N. R. Barlow/ W. Bhimji/ N. Chevaher/ P. J. Clark/ W. N. Cottmgham/ N. Dyce/ B. Foster/ 
Q ! C. Mackay/ D. Wallom/ F. F. Wilson/ K. Abe/ C. Hearty/ T. S. Mattison/ J. A. McKenna/ D. Thiessen/ 
Q ; S. JolIy/° A. K. McKemey/° V. E. Blinov/^ A. D. Bukin/^ D. A. Bukin/^ A. R. Buzykaev/^ V. B. Golubev/^ 
■ ^- ^- Ivanchenko/^ A. A. Korol/^ E. A. Kravchenko/^ A. P. Onuchin/^ A. A. Salnikov/^ S. I. Serednyakov/^ 
^ ; Yu. I. Skovpen/i V. I. Telnov/i A. N. Yushkov/^ D. Best/^ M. Chao/^ A. J. Lankford/^ M. Mandelkern/^ 
S. McMahon/2 D. P. Stoker/^ K. Arisaka/^ C. Buchanan/^ S. Chun/^ D. B. MacFarlane/-* S. PreU/^ 
Sh. Rahatlou/-* G. Raven/'' V. Sharma/'* C. Campagnari/^ B. Dahmes/^ P. A. Hart/^ N. Kuznetsova/^ 
S. L. Levy/5 Long/^ A. Lu/^ J. D. Richman/^ W. Verkerke/^ M. Witherell/^ S. Yellin/^ J. Beringer/^ 
>^ ■ D. E. Dorfan/6 A. M. Eisner/*^ A. A. Grillo/^ M. Grothe/^ C. A. Heusch/^ R. P. Johnson/^ W. S. Lockman/^ 
O ; T. Pulliani/6 H. Sadrozinski/^ T. Schalk/^ R. E. Schmitz/^ B. A. Schunim/^ A. Seiden/^ M. Turri/^ 

W. Walkowiak/6 D. C. WiUiams/^ M. G. Wilson/^ E. Chen/^ G. P. Dubois-Felsmann/^ A. Dvoretskii/^ 
D. G. Hitlin/^ S. Metzler/^ J. Oyang/^ F. C. Porter/^ A. Ryd/^ A. Samuel/^ M. Weaver/^ S. Yang/^ 
D ■ R. Y. Zhu/^ S. Devmal/8 T. L. Geld/^ S. Jayatilleke/^ G. Mancinelh/^ B. T. Meadows/^ M. D. SokolofF/^ 

^ T. Barillari/9 P. Bloom/^ M. O. Dima/^ S. Fahey/^ W. T. Ford/^ D. R. Johnson/^ U. Nauenberg/^ A. Olivas/^ 
D . P. Rankin/9 Roy/9 S. Sen/^ J. G. Smith/^ W. C. van Hoek/^ D. L. Wagner/^ J. Blouw/" J. L. Harton/° 
P^; M. Krishnamurthy/'^ A. Soffer/o W. H. Toki/o R. J. Wilson/" J. Zhang/" R. Aleksan/^ A. de Lesquen/^ 

1^ '■ S. Emery/1 A. Gaidot/^ S. F. Ganzhur/^ P.-F. Giraud/^ G. Hamel de Monchenault/^ W. Kozanecki/^ 

M. Langer/i G. W. London/^ B. Mayer/^ B. Serfass/^ G. Vasseur/^ Ch. Yeche/^ M. Zito/^ T. Brandt/^ 
J. Brose/2 T. Colberg/^ M. Dickopp/^ R. S. Dubitzky/^ A. Hauke/^ E. Maly/^ R. MiiUer-Pfefferkorn/^ 
^ i S. Otto/2 K. R. Schubert/2 R. Schwierz/^ B. Spaan/^ L. Wilden/^ D. Bernard/^ G. R. Bonneaud/^ 

5j ; F. Brochard/^ J. Cohen- Tanugi/^ S. Ferrag/^ E. Roussot/^ S. T'Jampens/^ Ch. Thiebaux/^ G. Vasileiadis/^ 
M. Verderi/3 A. Anjomshoaa/-* R. Bernet/* A. Khan/^ D. Lavin/"* F. Muheim/" S. Playfer/^ J. E. Swain/" 
J. Tinslay/"* M. Falbo/^ C. Borean/^ C. Bozzi/*' S. Dittongo/^ L. Piemontese/^ E. Treadwell/"^ F. AnuUi/**'! 
R. Baldini-Ferroh/* A. Calcaterra/^ R. de Sangro/^ D. Falciai/^ G. Finocchiaro/* P. Patteri/® L M. Peruzzi/*'0 
M. Piccolo/* Y. Xie/8 A. Zallo/* S. Bagnasco/^ A. Buzzo/^ R. Contri/^ G. Crosetti/^ M. Lo Vetere/^ 
M. Macri/9 M. R. Monge/^ S. Passaggio/^ F. C. Pastore/^ C. Patrignani/^ M. G. Pia/^ E. Robutti/^ 
A. Santroni/9 S. Tosi/^ M. Morii/o R. Bartoldus/^ R. Hamihon/i U. Mallik/i J. Cochran/^ H. B. Crawley/^ 
P.-A. Fischer/2 J. Lamsa/^ W. T. Meyer/^ E. I. Rosenberg/^ G. Grosdidier/^ C. Hast/^ A. Hocker/^ 
H. M. Lacker/3 S. Laplace/^ V. Lepeltier/^ A. M. Lutz/^ S. Plaszczynski/^ M. H. Schune/^ S. Trincaz-Duvoid/^ 
G. Wormser/3 R. M. Bionta/"* V. Brigljevic/'' D. J. Lange/'* M. Mugge/'* K. van Bibber/-* D. M. Wright/-* 
M. CarroU/5 J. R. Fry/^ E. Gabathuler/^ R. Gamet/^ M. George/^ M. Kay/^ D. J. Payne/^ R. J. Sloane/^ 

C. Touramanis/5 M. L. AspinwaU/'^ D. A. Bowerman/'^ P. D. Dauncey/'' U. Egede/'^ L Eschrich/*^ 
N. J. W. Gunawardane/6 J. A. Nash/^ P. Sanders/^ D. Smith/*^ D. E. Azzopardi/^ J. J. Back/"^ P. Dixon/"^ 
P. F. Harrison/^ R. J. L. Potter/^ H. W. Shorthouse/^ P. Strothcr/^ P. B. Vidal/^ M. L Wilhams/^ G. Cowan/* 
S. George/* M. G. Green/* A. Kurup/* C. E. Marker/* P. McGrath/* T. R. McMahon/* S. Ricciardi/* 
F. Salvatore/* L Scott/* G. Vaitsas/* D. Brown/^ C. L. Davis/^ J. AUison/o R. J. Barlow /o J. T. Boyd/° 
A. C. Forti/° J. Fullwood/° F. Jackson/o G. D. Lafferty/" N. Savvas/" E. T. Simopoulos/° J. H. Weatherall/" 



2 



A. Farbin,4i A. Jawahery ,"1 V. LiUard,''! J. Olsen,"! D. A. Roberts ,"1 J. R. Schieck G. Blaylock ,"2 

C. Dallapiccola,'*^ K. T. Flood,'*^ S. S. Hertzbach,42 R. Kofler,42 V. G. Koptchev ,^2 T. B. Moore,'*^ H. Staengle,'*^ 
S. Willocq,''2 B. Brau,''3 R. Cowan ^3 G. Sciolla,'*^ F. Taylor ,^3 R. K. Yamamoto,43 jyj^ Milek,^^ p. M. Patel,'*'* 

F. Palombo,'^^ J. M. Bauer,'*^ L. Cremaldi,'*^ V. Eschenburg,'*^ Kroeger,'*^ j Reidy,-*^ D. A. Sanders,"*^ 

D. J. Summers,'*'^ J. P. Martin,^^ J. Y. Nief,^^ R. Seitz,^^ p. Taras,"*^ V. Zacek H. Nicholson,48 C. S. Sutton,-** 
C. Cartaro,49 N. Cavallo,49^[| G. De Nardo,*^ F. Fabozzi,49 C. Gatto,''^ L. Lista,''^ P. Paolucci,^^ D. Piccolo,49 

C. Sciacca,^^ J. M. LoSecco,^" J. R. G. Alsmiller,^i T. A. Gabriel,^* T. Handler,^* J. Brau,^^ j^^ Fi-ey,^^ 
M. Iwasaki,52 ^ Sinev,^^ D. Strom,^^ Colecchia,53 F. Dal Corso,53 a. Dorigo,53 F. Galeazzi,53 Margoni,^3 

G. Michelon,53 jyj^ Morandin,53 M. Posocco,^3 Rotondo,^3 Simonetto,53 R. Stroili,53 E. Torassa,^3 

C. Voci,53 M. Benayoun,54 H. Briand,'^'' J. Chauveau,^'' P. David,^'* Ch. de la Vaissiere,^'* L. Del Buono,54 
O. Hamon,54 F. Le Diberder,^* Ph. Leruste,^^ J. Ocariz,54 L. Roos,^^ J. Stark,54 S. Versille,^'* P. F. Manfredi,^^ 
V. Re,55 V. Speziali,55 E. D. Frank,^^ l. Gladney,^^ Q. H. Guo,^^ J. Panetta,^^ C. Angelini," G. Batignani,^^ 

S. Bettarini,^'^ M. Bondioli," M. Carpinelli," F. Forti,^^ M. A. Giorgi,^^ A. Lusiani,^'^ F. Martinez- Vidal,^^ 
M. Morganti,^^ N. Neri,'^^ E. Paoloni,^'^ M. Rama,'" G. Rizzo," F. Sandrelli," G. Simi," G. Triggiani,^^ 
J. Walsli,57 M. Haire,58 D. Judd,^* K. Paick,^* L. Turnbull,^* D. E. Wagoner,^* J. Albert,^^ P. Elmer,^^ C. Lu,^^ 
K. T. McDonald,59 V. Miftakov,^^ S. F. SchafFner,^^ A. J. S. Smitli,^^ A. Tumanov,^^ E. W. Varnes,^^ G. Cavoto,60 

D. del Re,s° R. Faccini,!'^'^" F. Ferrarotto,^^ F. Ferroni,^^ E. Lamanna,60 E. Leonardi,'^'' M. A. Mazzoni,^" 

S. Morganti,60 G. Piredda,60 F. Safai Tehrani,^" M. Serra,60 C. Voena,s° S. Christ,^! R. Waldi,^! T. Adye,62 N. De 

Groot,®'^2 B_ Franek,62 N. I. Geddes,62 G. P. Gopal,62 g Xella,62 N. Copty,63 y Pui.ohit,63 H. Singli,'^3 
F. X. Yumiceva,^^ I. Adam,^^ P. L. Anthony,^'' D. Aston,*^" K. Baird,^^ N. Berger,^'* E. Bloom,^'* A. M. Boyarski,^'' 

F. Bulos,^'' G. Calderini,^'* M. R. Convery,^'^ D. P. Coupal,^'' D. H. Coward,^"* J. Dorfan,^^ W. Dunwoodie,^'' 
R. C. Field,^* T. Glanzman,64 G. L. Godfrey,64 Gowdy,^'' P. Grosso,^^ T. Haas,^'' T. Himel,^^ T. Hryn'ova,64 

M. E. Huffer,64 W. R. Innes,^'^ C. P. Jessop,^'* M. H. Kelsey,^'' P. Kim,^^ M. L. Kocian,64 U. Langenegger,^* 
D. W. G. S. Leith,64 S. Luitz,64 V. Luth,^* H. L. Lynch,64 H. Marsiske,^^ S. Menke,^'* R. Messner,^^ K. C. Moffeit,^^ 
R. Mount,64 D. R. Muller,^^ C. P. O'Grady,^'* V. E. Ozcan,^^ M. Pcrl,64 S. Petrak,^^ H. Quinn,*^" 

B. N. RatclifF,^'' S. H. Robertson,64 L. S. Rochester,*^" A. Roodman,64 T. Scliietinger,^^ R. H. Schindler,^^ 

J. Scliwiening,64 V. V. Serbo,^'^ A. Snyder,^* A. Soha,^'^ S. M. Spanier,64 J. Stelzer,^^ D. Su,^'' M. K. Sullivan,^^ 

H. A. Tanaka,^'' J. Va'vra,^'^ S. R. Wagner,^^ A. J. R. Weinstein,64 W. J. Wisniewski,^^ D. H. Wright,^^ 
C. C. Young,64 p. R. Burchat,^^ C. H. Cheng,^^ D. Kirkby,^^ T. 1. Meyer,^^ C. Roat,^^ R. Henderson,'^^' 

W. Bugg,67 H. Colin,67 A. W. Weidemann,^^ J. M. Izen,^^ I. Kitayama,^^ X. C. Lou,^* F. Bianchi,69 M. Bona,^^ 

D. Gamba,69 A. Smol,^^ L. Bosisio,™ G. Delia Ricca,™ L. Lanceri,™ P. Poropat,™ G. Vuagnin,^" R. S. Panvini,''! 

C. M. Brown,'^2 p,^ Jackson,^2 Kowalewski,'^2 Roney,^2 R. Band,'^3 p^ Charles, ^^3 g. Dasu,'^3 

F. Di Lodovico,'^3 a. M. Eichenbaum,'^3 jj 73 j ^ Johnson,'^3 pj^ 73 y. Pan,'^3 ^ Prepost,'^3 j j Scott,'^3 
S. J. Sekula,'^3 ^l. von Wimmersperg-Toeller,'^3 p^ -^^^3 2. Yu,^3 M. B. Kordich,74 and H. NeaF'* 

(The BABAR Collaboration) 

'Laboratoire de Physique des Particules, F-74941 Annecy-le-Vieux, France 
^Universita di Bari, Dipartimento di Fisica and INFN, F70126 Bari, Italy 
'''institute of High Energy Physics, Beijing 100039, China 
''' University of Bergen, Inst, of Physics, N-5007 Bergen, Norway 
^Lawrence Berkeley National Laboratory and University of California, Berkeley, CA 94720, USA 
^University of Birmingham, Birmingham, B15 2TT, United Kingdom 
''Ruhr Universitdt Bochum, Institut fiir Experimentalphysik 1, D-44^80 Bochum, Cermany 
'^University of Bristol, Bristol BS8 ITL, United Kingdom 
^University of British Columbia, Vancouver, BC, Canada V6T IZl 
^''Brunei University, Uxbridge, Middlesex UBS 3PII, United Kingdom 
Budker Institute of Nuclear Physics, Novosibirsk 630090, Russia 
University of California at Irvine, Irvine, CA 92697, USA 
''''University of California at Los Angeles, Los Angeles, CA 90024, USA 
'^ University of California at San Diego, La Jolla, CA 92093, USA 
"'University of California at Santa Barbara, Santa Barbara, CA 93106, USA 
'"University of California at Santa Cruz, Institute for Particle Physics, Santa Cruz, CA 95064, USA 
'"'California Institute of Technology, Pasadena, CA 91125, USA 
"^University of Cincinnati, Cincinnati, OH 45221, USA 
"'University of Colorado, Boulder, CO 80309, USA 
'"^^ Colorado State University, Fort Collins, CO 80523, USA 



3 



DAPNIA, Commissariat a I'Energie Atomique/Saclay, F-91191 Gif-sur-Yvette, France 
'"^^ Technische Universitdt Dresden, Institut fur Kern- und Teilchenphystk, D-01062, Dresden, Germany 
'^'^Ecole Polytechnique, F-91128 Palaiseau, France 
University of Edinburgh, Edinburgh EH9 3JZ, United Kingdom 
^^Elon University, Elon University, NC 27244-2010, USA 
Universitd di Ferrara, Dipartimento di Fisica and INFN, 1-44 100 Ferrara, Italy 
^''Florida A&M University, Tallahassee, FL 32307, USA 
^'^ Laboratori Nazionali di Frascati dell'INFN, 1-00044 Frascati, Italy 
^^Universitd di Geneva, Dipartimento di Fisica and INFN, I-I6I46 Genova, Italy 
^"Harvard University, Cambridge, MA 02138, USA 
University of Iowa, Iowa City, lA 52242, USA 
^^lowa State University, Ames, I A 50011-3160, USA 
Laboratoire de I'Accelerateur Lineaire, F-91898 Orsay, France 
'^^ Lawrence Livermore National Laboratory, Livermore, CA 94550, USA 
®* University of Liverpool, Liverpool L69 3BX, United Kingdom 
University of London, Imperial College, London, SW7 2BW, United Kingdom 
^^Queen Mary, University of London, El 4^^, United Kingdom 
* University of London, Royal Holloway and Bedford New College, Egham, Surrey TW20 OEX, United Kingdom 

University of Louisville, Louisville, KY 40292, USA 
^'^ University of Manchester, Manchester MIS 9PL, United Kingdom 
^' University of Maryland, College Park, MD 20742, USA 
University of Massachusetts, Amherst, MA 01003, USA 
^^Massachusetts Institute of Technology, Laboratory for Nuclear Science, Cambridge, MA 02139, USA 
^^McGill University, Montreal, QC, Canada H3A 2T8 
^^Universitd di Milano, Dipartimento di Fisica and INFN, 1-20133 Milano, Italy 
University of Mississippi, University, MS 38677, USA 
'^^ Universite de Montreal, Laboratoire Rene J. A. Levesque, Montreal, QC, Canada H3C 3J7 
Mount Holyoke College, South Hadley, MA 01075, USA 
^^Universitd di Napoli Federico II, Dipartimento di Scienze Fisiche and INFN, 1-80126, Napoli, Italy 
^° University of Notre Dame, Notre Dame, IN 46556, USA 
"Oak Ridge National Laboratory, Oak Ridge, TN 37831, USA 
■'''•^University of Oregon, Eugene, OR 97403, USA 
'^'^ Universitd di Padova, Dipartimento di Fisica and INFN, 1-35131 Padova, Italy 
Universites Paris VI et VII, Lab de Physique Nucleaire H. E., F-75252 Paris, France 
^^Universita di Pavia, Dipartimento di Elettronica and INFN, 1-27100 Pavia, Italy 
University of Pennsylvania, Philadelphia, PA 19104, USA 
^"^ Universita di Pisa, Scuola Normale Superiore and INFN, 1-56010 Pisa, Italy 
^'^Pravne View A&M University, Prairie View, TX 77446, USA 
'^'"Princeton University, Princeton, NJ 08544; USA 
''" Universitd di Roma La Sapienza, Dipartimento di Fisica and INFN, 1-00185 Roma, Italy 
Universitdt Rostock, D-18051 Rostock, Germany 
Rutherford Appleton Laboratory, Chilton, Didcot, Oxon, 0X11 OQX, United Kingdom 
University of South Carolina, Columbia, SC 29208, USA 
^^Stanford Linear Accelerator Center, Stanford, CA 94309, USA 
Stanford University, Stanford, CA 94305-4060, USA 
^^TRIUMF, Vancouver, BC, Canada V6T 2A3 
^''University of Tennessee, Knoxville, TN 37996, USA 
University of Texas at Dallas, Richardson, TX 75083, USA 
Universitd di Torino, Dipartimento di Fiscia Sperimentale and INFN, 1-10125 Torino, Italy 
'"Universitd di Trieste, Dipartimento di Fisica and INFN, 1-34127 Trieste, Italy 
''' Vanderbilt University, Nashville, TN 37235, USA 
University of Victoria, Victoria, BC, Canada V8W 3P6 
''^University of Wisconsin, Madison, WI 53706, USA 
''^Yale University, New Haven, CT 06511, USA 
(Dated: October 25, 2001) 

We present a measurement of the time-dependent CP-violating asymmetries in neutral B decays 
to the -K^-K^ CP cigcnstatc, and an updated measurement of the charge asymmetry in K^n^ 
decays. In a sample of 33 million T{4S) —^ BB decays collected with the BABAR detector at 
the SLAG PEP-II asymmetric B Factory, we find tt+tt" and 217 ± 18 K+n~ candidates 

and measure the asymmetry parameters /S^tt = O-OSlQ^g ± 0.11, Cnn = — 0.251q;47 ^ 0.14, and 
Aktt = —0.07 ± 0.08 ± 0.02, where the first error is statistical and the second is systematic. 
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In the Standard Model, all CP-violating effects arise 
from a single complex phase in the three-generation 
Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing ma- 
trix 1^ . One of the central questions in particle physics is 
whether this mechanism is sufficient to explain the pat- 
tern of CP violation observed in nature. Recent mea- 
surements of the parameter sin2/3 by the BABAR Q and 
BELLE 1^ Collaborations establish that CP symmetry 
is violated in the neutral i?-meson system. In addition to 
measuring sin2/3 more precisely, one of the primary goals 
of the i3-Factory experiments in the future will be to 
measure the remaining angles (a and 7) and sides of the 
Unitarity Triangle in order to further test whether the 
Standard Model description of CP violation is correct. 

The study of B decays to charmless hadronic two-body 
final states will play an increasingly important role in our 
understanding of CP violation. In the Standard Model, 
the time-dependent CP-violating asymmetry in the re- 
action B^ 7r+7r~ is related to the angle a. In addi- 
tion, observation of a significant rate asymmetry between 
B^ K^TT^ and 5° K~tt^ decays would be evidence 
for direct CP violation, and ratios of branching fractions 
for various tttt and Kir decay modes are sensitive to the 
angle 7. Finally, branching fraction measurements pro- 
vide critical tests of theoretical models that are needed 
to extract reliable information on CP violation from the 
experimental observables. 

The BABAR Collaboration recently reported measure- 
ments of branching fractions and charge asymmetries for 
several charmless two-body B decays using a data set of 
23 million BB pairs Q. In this paper, using a data sam- 
ple of approximately 33 million BB pairs, we report a 
measurement of time-dependent CP-violating asymme- 
tries in neutral B decays to the tt+tt" CP eigenstate 
and an updated measurement of the charge asymmetry 
in B^ K^TT~ decays. 

The time-dependent CP-violating asymmetry in the 
decay B'^ —> tt+tt" arises from interference between mix- 
ing and decay amplitudes, and interference between the 
tree and penguin decay amplitudes. A B^B^ pair pro- 
duced in T{AS) decay evolves in time in a coherent P- 
wave state until one of the two mesons decays. We re- 
construct a sample of B mesons {Bhh) decaying to the 
h'^h'~ final state, where h is a pion or kaon, and ex- 
amine the remaining charged particles in each event to 
"tag" the flavor of the other B meson (Stag)- The de- 
cay rate distribution /+ (/_) when h'^h'~ = tt+tt^ and 
Stag = B° (5°) is given by § 



and Ctttt are defined as 



/±(Ai) = 



\At\/T 



4t 



-[1 ± S'tttt sin(ATO(jAi) 

T C^^cos(AmdAi)], (1) 



where r is the i?" lifetime, Am^ is the B^B^ mixing 
frequency, and At — thh ~ ^tag is the time between the 
Bhh and -Stag decays. The CP-violating parameters S-m^ 



S-iTTT — ~ and Ctttt — 



i + iAr 



i + iAr 



(2) 



If the decay proceeds purely through the tree process 
b — > uW~ , the complex parameter A is directly related 
to CKM matrix elements, 



X{B 



v:hVtd\ fv:^Vub 



VudVZ 



(3) 



where we are assuming equal widths (AFs = 0) for the 
heavy and light mass eigenstates. Thus, at tree level in 
the Standard Model, |A| = 1 and ImX ~ sin2Q;, where 
a^s.vg[-VuV:jV^dV:,]. 

Recent theoretical estimates indicate that the contri- 
bution from the gluonic penguin amplitude can be sig- 
nificant §, 0, §. The process b dg carries the weak 
phase arg(V(^V(h), which can modify both the magni- 
tude and phase of A. Thus, in general, |A| 7^ 1 and 
2mX — |A|sin2acff, where acff depends on the magni- 
tudes and strong phases of the tree and penguin ampli- 
tudes. Several approaches have been proposed to obtain 
information on a in the presence of penguins ^, |^ . 

In this analysis, we extract signal and background 
yields for tt+tt^, K'^t:^, and K^K~ decays (lO|, and 
the amplitudes of the tttt sine (S'^tt) and cosine (C^tt) 
oscillation terms simultaneously from an unbinned maxi- 
mum likelihood fit. We parameterize the Kit component 
in terms of the total yield and the CP-violating charge 
asymmetry 



K- 



N 



K + 7 



K- 



(4) 



The data sample used in this analysis consists of 
33.7 fb~^ collected with the BABAR detector at the SLAC 
PEP-II storage ring between October 1999 and June 
2001. The PEP-II facility operates nominally at the 
T{AS) resonance, providing coUisions of 9.0 GeV electrons 
on 3.1 GeV positrons. The data set includes 30.4fb~^ col- 
lected in this configuration (on-resonance) and 3.3 fb~^ 
collected below the BB threshold (off-resonance) that are 
used for continuum background studies. 

A detailed description of the BABAR detector is pre- 
sented in Ref. jl|]. Charged particle (track) momenta 
are measured in a tracking system consisting of a 5-layer 
double-sided silicon vertex tracker (SVT) and a 40-layer 
drift chamber (DCH) filled with a gas mixture of he- 
lium and isobutane. The SVT and DCH operate within 
a 1.5T superconducting solenoidal magnet. The typical 
decay vertex resolution for fully reconstructed B decays 
is approximately 65 /im along the center-of-mass (CM) 
boost direction. Photons are detected in an electromag- 
netic calorimeter (EMC) consisting of 6580 CsI(Tl) crys- 
tals arranged in barrel and forward endcap subdetectors. 
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The flux return for the solenoid is composed of multiple 
layers of iron and resistive plate chambers for the identi- 
fication of muons and long-lived neutral hadrons. 

Tracks from the Bhh decay are identified as pions or 
kaons by the Cherenkov angle 9c measured with a detec- 
tor of internally reflected Cherenkov light (DIRC). The 
typical separation between pions and kaons varies from 
8(7 at 2 GeV/c to 2.5(T at 4GeV/c, where a is the aver- 
age resolution on 9c. Lower momentum kaons used in 
B flavor tagging are identified with a combination of 9c 
(for momenta down to 0.7GeV/c) and measurements of 
ionization energy loss dE/dx in the DCH and SVT. 

Hadronic events are selected based on track multiplic- 
ity and event topology. We require at least three tracks 
in the laboratory polar angle region 0.41 < 6'iab < 2.54 
satisfying the following requirements: transverse momen- 
tum greater than lOOMeV/c, at least 12 DCH hits, and 
originating from the interaction point within 10 cm in z 
and 1.5cm in i — Lp ]l^ . Residual two-prong events from 
the reaction e+e" {I = e,ii,T) are suppressed by 

requiring the ratio of Fox- Wolfram moments H2 / Hq ||1^ 
to be less than 0.95 and the sphericity jlj] of the event 
to be greater than 0.01. 

Candidate Bhh decays are reconstructed from pairs of 
oppositely-charged tracks forming a good quality vertex, 
where the Bhh four- vector is calculated assuming the pion 
mass for both tracks. We require each track to have 
an associated 9c measurement with a minimum of six 
Cherenkov photons above background, where the average 
is approximately 30 for both pions and kaons. Protons 
are rejected based on 9c and electrons are rejected based 
on dE/dx, shower shape in the EMC, and the ratio of 
shower energy and track momentum. Background from 
the reaction e+e~ ^ qq (q — u, d, s, c) is suppressed 
by removing jet-like events from the sample: we define 
the CM angle 9s between the sphericity axes of the B 
candidate and the remaining tracks and photons in the 
event, and require |cos0s| < 0.8, which removes 83% of 
the background. The total efficiency on signal events for 
all of the above selection is approximately 38%. 

We define a beam-energy substituted mass tties = 
\/E'^ — 'Pb- '^1^'^ candidate energy is defined as i?b — 
(s/2 -\- -pi ■ -pb)/ Ei, where ^/s and Ei are the total ener- 
gies of the e+e~ system in the CM and laboratory frames, 
respectively, and pi and ps are the momentum vectors 
in the laboratory frame of the e^e~ system and the Bhh 
candidate, respectively. Signal events are Gaussian dis- 
tributed in mEs with a mean near the B mass and a 
resolution of 2.6MeV/c^, dominated by the beam energy 
spread. The background shape is parameterized by a 
threshold function |l^ with a fixed endpoint given by 
the average beam energy. 

We define a second kinematic variable Ai? as the dif- 
ference between the energy of the Bhh candidate in the 
CM frame and \/s/2. The IS.E distribution is peaked 
near zero for tt+tt" decays. For decays with one (two) 



kaons, the distribution is shifted relative to tttt on av- 
erage by — 45MeV (— 91MeV), respectively, where the 
exact separation depends on the laboratory momentum 
of the kaon(s). The resolution on Ai? for signal decays 
is approximately 26 MeV. The background is parameter- 
ized by a quadratic function. 

Candidate h'^h'~ pairs selected in the region 5.2 < 
mES < 5.3GeV/c^ and |A£;| < 0.15 GeV are used to ex- 
tract yields and CP-violating asymmetries with an un- 
binned maximum likelihood fit. The total number of 
events in the fit region satisfying all of the above criteria 
is 9741. A sideband region, defined as 5.20 < rriEs < 
5.26 GeV/c^ and | Ai;| < 0.42 GeV, is used to extract var- 
ious background parameters. 

The analysis method combines the techniques used to 
measure charmless two-body branching fractions Q and 
sin2/3 The primary issues in this analysis are deter- 
mination of the i?tag flavor, measurement of the distance 
A2 between the Bhh and Stag decay vertices, discrimi- 
nation of signal from background, identification of pions 
and kaons, and extraction of yields and CP asymmetries. 

To determine the flavor of the Stag meson we use the 
same _B-tagging algorithm used in the sin2/3 and B^-B^ 
mixing analyses. The algorithm relies on the corre- 
lation between the flavor of the h quark and the charge 
of the remaining tracks in the event after removal of the 
Bhh candidate. We define five mutually exclusive tag- 
ging categories: Lepton, Kaon, NTl, NT2, and Untagged. 
Lepton tags rely on primary electrons and muons from 
semileptonic B decays, while Kaon tags exploit the cor- 
relation in the process 6 — > c — *■ s between the net kaon 
charge and the charge of the b quark. The NTl and NT2 
categories are derived from a neural network that is sen- 
sitive to charge correlations between the parent B and 
unidentified leptons and kaons, soft pions, or the charge 
and momentum of the track with the highest CM mo- 
mentum. The addition of Untagged events improves the 
signal yield estimates and provides a larger sample for de- 
termining background shape parameters directly in the 
maximum likelihood fit. 

The quality of tagging is expressed in terms of the ef- 
fective efficiency Q — ^iD^, where is the fraction of 
events tagged in category i and the dilution Di = 1 — 2wi 
is related to the mistag fraction Wi. The statistical er- 
rors on Stttt and Ctttt are proportional to 1/^/Q. Ta- 
ble I summarizes the tagging performance in a data sam- 
ple i?flav of fully reconstructed neutral B decays into 
D(*^-h+ {h+ = 7r+,p+,a+) and J/^K*° {K*° K+n-) 
fiavor eigenstates We use the same tagging efficiencies 
and dilutions for signal tttt, Ktt, and KK decays. Sepa- 
rate background tagging efficiencies for each species are 
obtained from a fit to the h~^h'~ on-resonance sideband 
data and reported in Table ||. 

The time difference At is obtained from the measured 
distance between the z position of the Bhh and Stag de- 
cay vertices and the known boost of the e'^e~ system. 
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TABLE I; Tagging efficiency e, average dilution D = 
1/2 {Dgo + Dgo), dilution difference AD — Dgo — Dga, and 
effective tagging efficiency Q for signal events in each tagging 
category. 

Category e{%) D {%) AD {%) Q {%) 

Lepton 11.0 ±0.3 82.3 ±2.7 -2.1 ±4.5 7.5 ± 0.5 

Kaon 35.8 ±0.5 64.8 ±2.0 3.5 ±3.1 15.0 ± 1.0 

NTl 8.0 ±0.3 55.6 ±4.2 -12.1 ± 6.7 2.5 ± 0.4 

NT2 13.9 ±0.4 30.2 ±3.8 9.0 ± 5.7 1.3 ± 0.3 
Untagged 31.3 ± 0.5 - 



Total Q 



26.3 ± 1.2 



TABLE IL Tagging efficiencies (%) for background events in 
each species. 

Category e(7r7r) e{KT7) e{KK) 

Lepton 1.0 ±0.1 1.0 ± 0.1 1.5 ± 0.2 

Kaon 26.0 ± 0.4 33.1 ± 0.6 23.5 ± 0.7 

NTl 6.6 ±0.2 5.4 ±0.3 6.9 ± 0.4 

NT2 17.6 ±0.4 15.3 ±0.5 19.7 ± 0.6 

Untagged 48.9 ± 0.7 45.2 ± 0.6 48.3 ± 0.8 



The z position of the -Stag vertex is determined with an 
iterative procedure that removes tracks with a large con- 
tribution to the total [@; • An additional constraint 
is constructed from the three-momentum and vertex po- 
sition of the Bfih candidate, and the average e"'~e~ in- 
teraction point and boost. The typical Az resolution is 
180 /zm. We require |At| < 17ps and 0.3 < cr/^t < 3. Ops, 
where a At is the error from the vertex fit. The resolution 
function for signal candidates is a sum of three Gaussians, 
identical to the one described in Ref. [||, with parame- 
ters determined from a fit to the Bfiav sample (including 
events in all five tagging categories). The background 
resolution function is parameterized as the sum of three 
Gaussians, with the parameters determined from a fit to 
the h~^h'~ on- resonance sideband data. 

The data sample used in the fit contains 97% back- 
ground, mostly due to random combinations of tracks 



produced in 



qq events. Discrimination of signal 



from background in the maximum likelihood fit is en- 
hanced by the use of a Fisher discriminant J- Q| . The dis- 
criminating variables are constructed from the scalar sum 
of the CM momenta of all tracks and photons (excluding 
tracks from the Bhh candidate) entering nine two-sided 
10-degree concentric cones centered on the thrust axis 
of the Bhh candidate. The distribution of for signal 
events is parameterized as a single Gaussian, with pa- 
rameters determined from Monte Carlo simulated decays 
and validated with B~ —> D'^Tr~ decays reconstructed in 



data. The background shape is parameterized as the sum 
of two Gaussians, with parameters determined directly in 
the maximum likelihood fit. 

Identification of h'^h'~ tracks as pions or kaons is ac- 
complished with the Cherenkov angle measurement from 
the DIRC. We construct Gaussian probability density 
functions (PDFs) from the difference between measured 
and expected values of 9c for the pion or kaon hypoth- 
esis, normalized by the resolution. The DIRC perfor- 
mance is parameterized using a sample of D*"*" D^Tr~^ , 
K~TT^ decays reconstructed in data. Within the 
statistical precision of the control sample (approximately 
10^ events), we find similar response for positively and 
negatively charged tracks and use a single parameteriza- 
tion for both. 

We use an unbinned extended maximum likelihood fit 
to extract yields and CP parameters from the B^h sam- 
ple. The likelihood for candidate j tagged in category 
c is obtained by summing the product of event yield ti^, 
tagging efficiency Ci^c, and probability Vi^c over the eight 
possible signal and background hypotheses i (referring to 
TTTT, K^TT~ , K~TT^, aud KK decays). 



£c = exp - ^ n^Ci^c Y]. 



E 



(5) 



For the K^ir^ hypotheses, the yield is parameterized as 
n-i = Nktt (1 ± Ak^) /2, where Nk^t = +Nk+-„- ■ 

We fix the tagging efficiencies to the values in Ta- 
bles I and H The probabilities Vi^c are evaluated as the 
product of PDFs for each of the independent variables 
Xj — {rriEs, ^E,J^,6'^ ,d~ , At}, where 9^ and 6~ are 
the Cherenkov angles for the positively and negatively 
charged tracks. The total likelihood £ is the product of 
likelihoods for each tagging category and the free param- 
eters are determined by minimizing the quantity In C. 

The At PDF for signal tt+tt^ decays is given by Eq. [|, 
modified to include the dilution and dilution difference 
for each tagging category, and convolved with the sig- 
nal resolution function. The At PDF for signal Ktt 
events takes into account B^-B^ mixing, depending on 
the charge of the kaon and the flavor of i^tag- We param- 
eterize i?" — > K^K~ decays as an exponential convolved 
with the resolution function. 

There are 18 free parameters in the fit. In addition 
to the CP- violating parameters S'tttt , C-mr , and Aktt , the 
fit determines signal and background yields (six param- 
eters), the background Kn charge asymmetry, and eight 
parameters describing the background shapes in rriEs, 
AE, and T. We fix r and Anid to the world-average 
values |l7| . 

In a sample of 33 million BB pairs we find 65l|[^ tttt, 
217 ± 18 Ktt, and 4.3l^ j KK events. These yields 
are consistent with the branching fractions reported in 
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TABLE III: Central values and 90% C.L. intervals for S^^, 
Ctttt, and Aktv from the maximum likelihood fit. 

Parameter Central Value 90% C.L. Interval 

0.03l2 j^ ± 0.11 [-0.89, +0.85] 
C^^ -0.25l°;*? ± 0.14 [-1.0, +0.47] 
Ak^ -0.07 + 0.08 + 0.02 [-0.21, +0.07] 



Ref. 1^, as well as measurements from other experi- 
ments |l^. The results for CP- violating asymme- 
tries are summarized in Table [II. Statistical errors cor- 
respond to unit change in ; 



-21n(£). For each pa- 
rameter, we also calculate the 90% confidence level (C.L.) 
interval corresponding to a change in of 2.69, and tak- 
ing into account the systematic error. The correlation be- 
tween S-jTTT and Ctttt is —21%, while Akh is uncorrelated 
with either 8^^-,^ or Ctttt- 

Figure [| shows distributions of toes and Ai? for 
events enhanced in signal decays based on likelihood ra- 
tios. We define TZs\g = Y^gnsTs/ Y^^niVi and TZk = 
nkPk/ X^s ^sVs, where (^.) indicates a sum over sig- 
nal (all) hypotheses, and Vk indicates the probability for 
signal hypothesis k. The probabilities include the PDFs 
for Oc, T , and toes (A-E-) when plotting Ai? (toes). The 
selection is defined by optimizing the signal significance 
with respect to 7?.sig and 'R-k- The solid curve in each 
plot represents the fit projection after correcting for the 
efficiency of the additional selection (approximately 55% 
for TTTT and 85% for K-k). 

Figure 1^ shows the Ai distributions and the asymmetry 
X.(Ai) = (iVBo(At)-iV5o(Ai))/(iVso(Ai) + Argo(A0) 
for tagged events enhanced in signal tttt decays. The se- 
lection procedure is the same as Fig. |l|, with the likeli- 
hoods defined including the PDFs for Qc-, wes, and 
Ai?. Approximately 24 tttt, 22 gg, and 5 K-k events sat- 
isfy the selection. 

Systematic uncertainties on S'^tt, Ct^t^, and Ak-k arise 
primarily from imperfect knowledge of the PDF shapes 
and uncertainties on tagging efficiencies, dilutions, r, and 
Atoc;. The total systematic error is calculated as the sum 
in quadrature of the individual uncertainties. The error 
on Ak-k is dominated by uncertainty in the mean of the 
A£' PDF (0.01) and possible charge bias in track and 6'c 
reconstruction (0.01) Errors on S-^: and C^^r are 

dominated by the parameterization of Af resolution for 
signal and background (« 0.07 for S'tttt, ~ 0.03 for C^tt), 
tagging (0.05), and, for Ctttt only, the mean of the Ai? 
PDF (0.1). 

Extensive studies were performed to validate the fit 
technique. A large ensemble of Monte Carlo pseudo- 
experiments was generated from the nominal PDFs with 
the statistics observed in the full data set. Parameter 
errors and the maximum value of the likelihood obtained 
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FIG. 1: Distributions of wes and A_E (unshaded histograms) 
for events enhanced in signal (a), (b) tttt and (c), (d) K-k 
decays based on the likelihood ratio selection described in 
the text. Solid curves represent projections of the maximum 
likelihood fit result after accounting for the efficiency of the 
additional selection, while dashed curves represent qq and 
TTTT ^ K-E cross-feed background. Shaded histograms show 
the subset of events that are tagged. 
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FIG. 2: Distributions of for events enhanced in signal tttt 
decays based on the likelihood ratio selection described in the 
text. Figures (a) and (b) show events (points with errors) 
with Stag = or _B". Solid curves represent projections 
of the maximum likelihood fit, dashed curves represent the 
sum of qq and K-e background events, and the shaded region 
represents the contribution from signal tttt events. Figure (c) 
shows A-n-n^l^t) for data (points with errors), as well as fit 
projections for signal and background events (solid curve), 
and signal events only (dashed curve). 



in the data fit are all consistent with expectations based 
on these pseudo-experiments, and all free parameters are 
unbiased. We have checked that consistent results are 
obtained when separating events by -Stag flavor. As a 



8 



validation of the At parameterization in data, we fit 
the full data set to simultaneously extract yields, back- 
ground parameters, t, Am^;, S.^t^, and C-^Tr- We find 
T = (1.52±0.12)ps and Am^ = (0.54± 0.09)?ips^i, and 
all other parameters are consistent with the nominal fit. 

In summary, we have presented a measurement of time- 
dependent CP-violating asymmetries in —> tt+tt" de- 
cays and an updated measurement of the charge asym- 
metry Aktt- The latter is consistent with our previous 
result reported in Rcf. Q|, as well as results from other 
experiments . We observe no evidence for direct 

CP violation in the Ktt mode and determine a 90% C.L. 
interval excluding a significant part of the allowed region. 
Although the current measurements of and Ctttt do 
not significantly constrain the Unitarity Triangle, with 
the addition of more data and further improvements in 
detector performance and analysis techniques, future re- 
sults will yield important information about CP violation 
in the i3-meson system. 
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